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ABSTRACT 

By combining complementary monitoring observations spanning long, medium and short time scales, 
we have constructed power spectral densities (PSDs) of six Seyfert 1 galaxies. These PSDs span >4 
£C) , orders of magnitude in temporal frequency, sampling variations on time scales ranging from tens of 

minutes to over a year. In at least four cases, the PSD shows a "break," a significant departure from a 
power law, typically on time scales of order a few days. This is similar to the behavior of Galactic X-ray 
binaries (XRBs), lower mass compact systems with breaks on time scales of seconds. NGC 3783 shows 
tentative evidence for a doubly-broken power law, a feature that until now has only been seen in the 
(much better-defined) PSDs of low-state XRBs. It is also interesting that (when one previously-observed 
object is added to make a small sample of seven), an apparently significant correlation is seen between 
the break time scale T and the putative black hole mass A/bh, while none is seen between break time 
scale and luminosity. The data are consistent with the linear relation T = Mbh/10 6 ' 5 M q ; extrapolation 
over 6-7 orders of magnitude is in reasonable agreement with XRBs. All of this strengthens the case for 
a physical similarity between Seyfert Is and XRBs. 

Subject headings: galaxies: active — galaxies: Seyfert — X-rays: galaxies 
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UMP02 to create a PSD survey of seven Seyfert 1 galaxies. 
The physical implications of this break is discussed and a 
comparison with the PSD and other variability properties 
of XRBs is made in § 5. A brief summary of these results 
is given in § 6. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Sampling 

The sampling pattern employed herein (based on EN99) 
covers complementary long, medium and short time scales 
with progressively shorter and denser relatively even sam- 
pling. This makes it possible to construct a PSD that 
covers the maximum temporal frequency range while min- 
imizing the amount of telescope time needed. It also re- 
quires the assumption that the variations are stationary or 
at most mildly non-stationary; this assumption is tested in 
§ 3.4. The specific sampling patterns are described below 
and tabulated in Table 2. 

• Long (months-year) time scales: For all six tar- 
gets, RXTE obtained nearly even sampling every 
~4.3 d (~64 orbits) over a period of >3 years. In 
NGC 5548, observations in 1996-1999 every -14 d 
were also included. Observations lasted only a frac- 
tion of an orbit (—1 ksec) in all cases. 

• Medium (day— weeks) time scales: For four of 
the objects, medium time scale sampling was ob- 
tained with RXTE by observing once every other 
orbit (3.2 hrs) for 256 epochs, spanning 34 days. 
For NGC 3783, this sampling was once every other 
orbit for 151 epochs (~20 days), and for NGC 3516, 
it was once every eight orbits for 136 days. Again, 
~1 ksec samples were used. 

• Short (hours) time scales: Quasi-continuous 
(—60-300 ksec) observations were obtained with 
XMM-Newton (for NGC 5548 and NGC 4151), 
Chandra (NGC 3783) and RXTE (NGC 3516). 
These observations were interrupted only by peri- 
ods of increased particle backgrounds (for RXTE 
and occasionally XMM-Newton) and were sensi- 
tive to >3% variations on time scales as short as 
—1 ksec. The other two sources (Akn 564 and 
Fairall 9) had only short (<35 ksec) uninterrupted 
XMM-Newton observations. 

2.2. RXTE data reduction 

All of the long- and medium-term data were obtained 
with RXTE using similar sampling patterns; data reduc- 
tion proceeded in a similar fashion for all these observa- 
tions, as well as the two NGC 3516 long-looks. These 
data were obtained with the proportional counter array 
(PCA), which consist of five identical collimated propor- 
tional counter units (PCUs; Swank 1998). For simplic- 
ity, data were collected only from those PCUs which did 
not suffer from repeated breakdown during on-source time 
(PCUs 0, 1, and 2 prior to 1998 December 23; PCUs and 
2 from 1998 December 23 until 2000 May 12; PCU 2 only 
after 2000 May 12). Count rates quoted in this paper are 
normalized to 1 PCU. Only PCA STAND ARD-2 data were 



considered. The data were reduced using standard extrac- 
tion methods and FTOOLS v4.2 software. Data were re- 
jected if they were gathered less than 10° from the Earth's 
limb, if they were obtained within 30 min after the satel- 
lite's passage through the SAA, if ELECTRON0 > 0.1 
(ELECTRON2 after 2000 May 12), or if the satellite's 
pointing offset was greater than 0?02. 

As the PCA has no simultaneous background monitor- 
ing capability, background data were estimated by using 
PCABACKEST v2.1b to generate model files based on the 
particle-induced background, SAA activity, and the diffuse 
X-ray background. This background subtraction is the 
dominant source of systematic error in RXTE AGN moni- 
toring data (e.g., EN99). Counts were extracted only from 
the topmost PCU layer to maximize the signal-to-noise ra- 
tio. All of the targets were faint (< 40 ct s _1 PCU -1 ), 
so the applicable 'L7-240' background models were used. 
Because the PCU gain settings changed three times since 
launch, the count rates were rescaled to a common gain 
epoch (gain epoch 3) by calibrating with several public 
archive Cas A observations. Light curves binned to 16 s 
were generated for all targets over the 2-10 keV bandpass, 
where the PCA is most sensitive and the systematic er- 
rors and background are best quantified. The data were 
then further binned as listed in Table 2, column 5; bins 
with less than 10 flux points were excluded from analysis. 
Standard errors were derived from the data in each orbital 
bin. Further details of RXTE data reduction can be found 
in EN99. 

2.3. XMM-Newton observations and data reduction 
2.3.1. NGC 5548 

XMM-Newton observed NGC 5548 for 97 ksec during 
revolution 290 from 2001 Jul 09 16:08:04 UT to 2001 Jul 10 
18:06:21 UT. Data from the high-throughput European 
Photon Imaging Camera (EPIC) instruments, the pn 
(Struder et al. 2001) and two MOS cameras (Turner et 
al. 2001) were used. The medium filters were used. Due 
to their higher count rate and signal-to-noise, only the 
pn data were used for PSD analysis; the MOS data are 
henceforth ignored. The pn camera was operated in Small 
Window mode; the readout time was 5.7 ms. 

A standard reduction of the raw pn data was done using 
the Science Analysis System. This involved the subtrac- 
tion of hot, dead, or flickering pixels, removal of events 
due to electronic noise and correction of event energies for 
charge transfer losses. The patterns used were 0-4 (single 
and double pixel events). The extraction region used was a 
circle of radius 40". The core was not excluded, as the level 
of pile-up was fairly small (<2%). The background count 
rate was extracted in 250 second bins using the same size 
region in the same chip, but off-source (>1')- Solar flare 
activity caused a tremendous increase in the soft proton 
flux during the final 13.5 ksec, rendering that section of 
data useless for PSD analysis. This resulted in 83.5 ksec 
of useful data; a background-subtracted light curve over 
the 2-10 keV band was initially extracted in 1 s bins. In 
a red-noise PSD derived from a continous light curve, the 
power amplitude at the highest temporal frequencies will 
be dominated by the constant power level contributed by 
Poisson noise, Pp sn , instead of the intrinsic source vari- 
ability. To ensure a high variability-to-noise PSD with the 



3 



intrinsic source varaibility power being greater than the 
Poisson noise power at all frequencies sampled, the data 
were binned at the time scale where the intrinsic variabil- 
ity was equal to Pp sn ; in the case of the NGC 5548 data; 
this time scale was 900 s. 

2.3.2. NGC 4151 

NGC 4151 was observed by XMM-Newton for 100 ksec 
during revolution 190 from 2000 Dec 22 02:48:22 UT to 
2000 Dec 23 05:29:59 UT. Due to two 8 ksec gaps which 
would have seriously complicated the PSD analysis, only 
the final 57 ksec was used. Again, the medium filter was 
used for all three EPIC instruments, which were operated 
in Full Frame mode. Data reduction over the 2-10 keV 
band proceeded in a similar manner to that of NGC 5548; 
again, due to its higher count rate and signal-to-noise, the 
pn data only were used for PSD analysis. The patterns 
used were 0-4 (single and double pixel events). The area 
of source extraction was a circle of radius 30"; background 
count rates were extracted over a region 4.4 times as large. 
Due to the low source flux, photon pile-up was not an issue. 
Background-subtracted light curves initially extracted in 
1 s bins were rebinned to 1100 s. 

2.4. Chandra observations and data reduction for 
NGC 3783 

Chandra observed NGC 3783 for 170 ksec each of five 
times during 2001 March-June using the High Energy 
Transmission Grating Spectrometer (HETGS; e.g., Mark- 
ert et al. 1994) with the Advanced CCD Imaging Spec- 
trometer (ACIS; e.g., Nousek et al. 1998; Plucinsky et 
al. 2002) as the detector. The data were reduced using 
CI AO software version 2.1.2 and its associated calibration 
data. Photons in the 2-10 keV band were extracted from 
the MEG and HEG dispersed first orders, as the zeroth 
order in the HETGS images of NGC 3783 suffers from 
heavy photon pile-up. The extraction width on the dis- 
persed spectrum was done with the CIAO default of 4.78". 
NGC 3783 is a point like source and no extended circum- 
nuclear emission was detected; the working assumption 
was that all photons in the dispersive arms were from the 
central source. The background count rate was negligible 
(several orders of magnitude smaller than the signal) in 
the dispersed arms, and was not subtracted. Further de- 
tails of the Chandra data reduction can be found in Kaspi 
et al. (2002). Light curves initially extracted in bins of 
32.41 s were truncated to equal lengths of 167.6 ksec and 
rebinned to 2000 s for a total of 84 data points in each 
light curve. 

2.5. The light curves 

The resulting light curves are shown in Figures 1-4. Fig- 
ure 1 shows the long- and medium-term data from RXTE 
monitoring, with the short-term NGC 3516 RXTE long- 
looks included and the locations of the XMM-Newton and 
Chandra long-looks marked. Figure 2 shows the long-term 
light curves resampled from the raw RXTE light curves. 
Figure 3 shows the resampled medium-term light curves, 
and Figure 4 shows the intensive long-looks. Table 2 sum- 
marizes the sampling parameters for each light curve, in- 
cluding the instrument, mean count rate, sampling inter- 
val, total points, and percent of missing data for each time 



scale for which there is even sampling suitable for PSD 
analysis. Also listed is F var , the square root of the ex- 
cess variance, a measure of the intrinsic variability am- 
plitude normalized to the mean as defined in, e.g., Edcl- 
son et al. (2002). We note that RXTE-PCA, Chandra- 
HETG, and XMM - Newton EPIC pn do not have pre- 
cisely identical responses over the 2-10 keV bandpass; 
the PCA response is somewhat harder (peaking closer to 
^5 keV) than that of the pn or Chandra-HETG (which 
peak closer to ~2 keV). It has been shown that measured 
high-frequency PSD slope is energy dependent in both 
XRBs and Seyferts (Nowak et al. 1999a, Nandra & Pa- 
padakis 2001), flattening as photon energy increases; our 
use of RXTE data at low frequencies and XMM-Newton or 
Chandra at high frequencies could thus potentially mimic 
a high-frequency break. However, we expect this effect 
to be relatively minor, as PSD slope typically flattens by 
only one or two tenths in power law slope for a doubling 
or tripling of photon energy (Nowak et al. 1999a, Nandra 
& Papadakis 2001), and XRB and PSD breaks typically 
involve a change in power law slope near 1. 

3. PSD CONSTRUCTION 

The PSD construction required several steps. First, 
PSDs were measured for each individual short, medium or 
long time scale light curve and then the individual PSDs 
were combined to form one high-dynamic range PSD for 
each target, as described in § 3.1. However, these com- 
bined PSDs have poorly determined errors (dominated by 
systematic effects such as aliasing and red-noise leak, as de- 
scribed in § 3.2). A Monte Carlo technique was employed 
to better estimate PSD fit parameters in the presence of 
these systematic errors and distortion effects, as described 
in § 3.3. 

3.1. Initial measurement 

First, light curves were linearly interpolated for missing 
data points, though such gaps were relatively rare. Light 
curves were not interpolated to a strictly even grid, since 
departures from ephemcris were almost always small (typ- 
ically only a few percent) and such an interpolation would 
have negligible impact on the resulting periodogram. Each 
light curve's mean was subtracted. 

The periodograms were constructed using a Discrete 
Fourier Transform (DFT; e.g., Oppenhcim & Shafer 1975; 
Brillingcr 1981). The power at each Fourier frequency 
/ = 1/D,2/D, ... 1/(2 AT) (the periodogram, where AT 
is the sampling time) was calculated using the fractional 
RMS squared normalization 

pu) = j^mt 

where D is the duration, N is the number of points, fi is the 
mean count rate of the data, and \F (f)\ 2 is the modulus- 
squared of the Fourier Transform of the light curve. The 
null, or rectangular, lag window was used; use of a tapered 
lag window for steep red-noise PSD data strongly increases 
the bias of the periodogram (see e.g., Priestley 1981). Fol- 
lowing Papadakis & Lawrence (1993a), the logarithm of 
the periodogram was binned every factor of 1.4 (0.15 in 
the logarithm) in temporal frequency, but the lowest fre- 
quency bins were widened to accommodate a minimum of 
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2 periodogram points. The constant level of power due to 
Poisson noise was not subtracted, but it was modelled in 
the Monte Carlo analysis. The individual PSDs were then 
combined to form broadband PSDs for each target. 

These high-quality broadband PSDs spanned an excep- 
tional range of time scales: four cover more than 4.4 
decades of temporal frequency and the other two cover 
3.6 decades. PSD measurement parameters, including the 
minimum and Nyquist frequencies sampled in each PSD 
segment, are given in Table 3. No renormalization in power 
amplitude of the individual PSDs was done. The resulting 
broadband PSDs are shown in Figure 5. Visual inspection 
reveals a variety of behavior; e.g., low- frequency flattening 
is readily evident in the PSDs of NGC 3783 and Akn 564. 

3.2. Aliasing and red-noise leak 

Although temporal analysis techniques such as the PSD 
are derived assuming continuous observations of infinite 
duration, such conditions cannot be attained in practice. 
These real-life Seyfert 1 sampling patterns S(t) span a 
limited duration D and have a shortest time resolution 
AT samp . As the Fourier transform of the observed data 
is a convolution of the Fourier transform of the underly- 
ing variability process with S(f), the 'sampling window,' 
power is transferred into the sampled frequency range 
from above /N yq = l/(2AT samp ) and below 1/D (red-noise 
leak), as explained below. 

If a red-noise PSD contains significant power below 
/min = 1/D, then there may be significant long-term 
trends present in the light curve. These trends can domi- 
nate the total variance of the observed light curve, whose 
measured PSD then contains additional power transferred 
from below / m i n . The power transferred by this red-noise 
leak process has a power-law slope which goes as f~ 2 (see 
details in van der Klis 1997). 

On short time scales, the discrete sampling causes alias- 
ing. For bins of width ATbi n evenly spaced AT samp apart 
(with ATbin <SC AT samp ), variations on time scales shorter 
than AT samp cannot be distinguished from longer-time 
scale variations, and power from above /Nyq is effectively 
added into the frequency range sampled by the PSD. 

Distortion effects arc present in, and hamper the in- 
terpretation of, the PSDs constructed above, but another 
serious problem with these PSDs is that the lowest tempo- 
ral frequency bins contain too few PSD points for normal 
errors to be assigned. By binning a sufficient number of 
periodogram estimates (>20 for logarithmically binned pe- 
riodograms; Papadakis & Lawrence 1993a), the averaged 
power will approach a Gaussian distribution, and the mean 
of the periodogram points in a frequency bin will tend to 
the power amplitude of the true underlying variability pro- 
cess, Ptrue(f)- 

3.3. Monte Carlo simulations 

To obtain a PSD shape with adequate errors, and to 
solve the problem of distortion effects, we use a version 
of Monte Carlo technique PSRESP first introduced by 
UMP02, based on a similar Monte Carlo technique by 
Done et al. (1992). This enables us to quantify the degree 
of low-frequency flattening and find a best-fit model PSD 
shape corresponding to the underlying variability process. 



The technique consists of simulating light curves from 
a given PSD model shape specified for testing, resampling 
these light curves and measuring their PSDs in the same 
manner as the observed light curves, forming an 'average 
model' broadband PSD which accounts for the distortion 
effects. Errors are assigned for all frequency bins, derived 
from the RMS spread of the individual simulated PSDs 
at a given frequency bin. Finally, the technique compares 
the model broadband PSD to that derived from the ob- 
served data using a fit statistic with the distribution esti- 
mated from the simulations to determine goodness of fit 
for the model. In this manner, a variety of underlying PSD 
model shapes can be tested against the data. We note, as a 
caveat, that the results of this technique are highly model 
dependent; specifically, the primary assumption govern- 
ing this entire process is that the broadband PSD model 
shape used for testing is an accurate representation of the 
underlying variability process Ptrue(f)- An outline of the 
Monte Carlo method is as follows (see UMP02 for further 
details): 

1. An underlying model PSD shape (a continuous func- 
tion) is specified for testing. The initial normalization is 
arbitrary. For each long-, medium-, and short-term PSD 
segment, the algorithm of Timmer & Konig (1995) is used 
to generate N light curves (where N is at least 100). The 
time resolution AT s ; m is chosen to be 0.1AT samp . To en- 
sure that the light curves account for variability on time 
scales longer than the observed light curve (contain a red- 
noise leak contribution if necessary), one very long segment 
of length N x D is simulated; this is then broken up into 
N light curves, each of duration D. 

2. The light curves are resampled in the same man- 
ner as the observed light curves as follows: For long- and 
medium-term simulated data, every tenth point is selected 
to degrade the resolution from AT s ; m = 0.1AT samp to 
AT samp ; this correctly accounts for most of the total alias- 
ing, specifically that portion of the aliased power which is 
due to variations on time scales from /Nyq to l/(2AT s i m ). 
For continuous long-look simulated data, the data are av- 
eraged over every ten points. The uncertainty in the quan- 
tity of aliased power in red-noise PSD is dominated by the 
uncertainty in the amount of power at frequencies just 
above /Nyq, which is why this resampling procedure must 
be used. Time bins with missing flux in the observed light 
curves are linearly interpolated. 

3. PSDs are constructed from the simulated light curves 
in the same manner as the observed PSDs, using the same 
normalization and frequency binning. For each segment, 
the model average PSD P s im(/) is calculated from the N 
individual PSDs; error bars A P s im(/) equal to the RMS 
spread of the individual PSDs at each frequency bin are 
assigned. 

4. While the resampling approach in step 2 accounts for 
much of the total aliasing in red-noise PSDs, there remains 
aliasing due to variations on time scales from AT s i m down 
to ATbin- This second, smaller quantity of aliased power 
is approximated by adding a constant level of power to 
all frequencies sampled in a given PSD. Since variations 
on time scales shorter than ~ 2ATbi n will likely not con- 
tribute to aliasing, this quantity of aliased power can be 
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approximated by the analytical expression 

1 r (2AT biD r 1 



alias 



/Nyq — /min J(2AT aim )- 1 



P(f)df 



(e.g., see UMP02 § 3.3) and added to the model average 
PSD and the N individual PSDs. The resulting PSDs 
then account fully for both red-noise leak and both alias- 
ing components. 

5. The goodness of fit is determined. Comparison of 
a model-average PSD to the observed PSD is not pos- 
sible in the "traditional" x 2 sense, since, as previously 
discussed, error bars assigned to the observed PSD are 
not strictly Gaussian. Instead, a new statistic Xdist i s de- 
fined to compare the average distorted model PSD to the 
observed data, using the well-determined errors from the 
model: 

2 _ V (P S im(/)-Pobs(/)) 2 

Xdist y (AP sim (/)) 2 

The total Xdist value is calculated by summing over all 
frequency bins in all PSD segments. For the targets with 
multiple short-term PSDs (NGC 3516 and NGC 3783), we 
average the values of short-term Xdist to avoid unnecessar- 
ily weighting the fit towards the highest frequencies. 

The best- fitting normalization of the P s i m (f) is found by 
renormalizing all segments of the P s i m (f) by the same fac- 
tor until the total Xdist value is minimized. The power con- 
tribution from Poisson noise is added to the model PSDs 
during this step; Pp sn is calculated via 2( / u + _B)/ y u 2 , where 
fj, and B are the total source and background count rates, 
respectively. For non-continuously observed light curves, 
this must be multiplied by AT samp /ATbin- The observed- 
Xdist va l ue is calculated. 

The goodness of fit is then determined as follows: 10,000 
combinations the of long-, medium-, and short-term sim- 
ulated PSDs arc randomly selected to model the Xdist dis- 
tribution for each given PSD model. For each combina- 
tion, the value of simulated-Xd; s t (comparing P s i m (/) to 
-Psim(/)) is calculated. These 10,000 measurements of the 
simulatcd-Xdi s t distribution arc sorted into ascending or- 
der. The probability that the model PSD can be rejected 
is given by the percentile of the 10,000 Xdist exceeded by 
the value of observed Xdist- 

6. The above steps are repeated to test a range of 
PSD model shapes by stepping through a range of high- 
frequency slopes and break frequencies. 

4. RESULTS 

The results of the Monte Carlo analysis are presented 
for several simple PSD model shapes to quantifying the 
degree of flattening, if any, towards low temporal fre- 
quencies. First, unbroken PSD models are tested (§ 4.1). 
Then, models incorporating a single PSD break are tested 
(§ 4.2). We do not test for quasi-periodic oscillations 
(QPO); though they are routinely seen in XRB PSDs, 
there is no obvious indication of QPO in the observed light 
curves or resultant PSDs, and there has been no convinc- 
ing evidence to date for deterministic behavior in AGN 
light curves or PSDs (e.g., Lawrence et al. 1987, McHardy 
& Czerny 1987) as previous claims of QPO in AGN (e.g., 
Papadakis & Lawrence 1993b) have all either been refuted 



or lack confirmation. Preliminary evidence for a doubly- 
broken power law model fit to the PSD of NGC 3783 is 
presented in § 4.3. Finally, a test of one of the assumptions 
governing the PSD construction, the assumption of PSD 
stationarity, is presented in § 4.4. 

4.1. Unbroken power law models 

The simplest model tested was an unbroken power law, 
of the form P(f) = j4 (/// ) - ^, where the normalization 
A a is the PSD amplitude at / = f , and (3 is the power law 
slope. (The constant level of power from Poisson noise is 
added to each simulated PSD, but as this value is different 
for each long-, medium-, and short-term PSD segment, it 
is not explicitly listed here.) For all targets, the model 
was tested by stepping through the range of (3 from 1.0 to 
4.0 in increments of 0.05; for Akn 564 and NGC 3783, the 
two targets with relatively flatter observed PSDs, values 
of (3 down to 0.0 were additionally tested. Five hundred 
simulated PSDs were used to determine the average model 
PSD, and 10 4 randomly selected sets of PSDs were used 
to probe the simulated Xdist distribution. 

The best-fitting models and simulated data are plotted 
in / x Pf space in Figure 6 for all targets. The best- 
fitting values of (3 and A a , along with the corresponding 
likelihood of acceptance (defined as 1 - R, where R is the 
rejection confidence) are summarized in Table 4. The er- 
rors on (3 correspond to values la above the likelihood of 
acceptance for the best-fit value on a Gaussian probability 
distribution (i.e., the amount (3 needs to change for the fit 
to be less likely by la). While we have not yet rigorously 
proven the validity of these errors, this method does give 
very reasonable values. The errors on A D are determined 
from the RMS spread of the 10 4 randomly selected sets 
of simulated PSDs. The low likelihood of acceptance for 
NGC 3783, for instance, indicates that an unbroken power 
law is an inaccurate description of the data. 

4.2. Singly-broken power law models 

To test for the presence of a PSD break, we employed a 
broken power law model of the form 

A{f/f c )~\ f<f c 
A(f/f c )- , f>f c 
where the normalization A is the PSD amplitude at the 
break frequency f c , (3 is the high frequency power law 
slope, and 7 is the low frequency power law slope, with 
the constraint 7 < (3. The range of slopes tested was 
(3 = 0.0 — 4.0 in increments of 0.05. Break frequencies 
were tested in the log from -8.0 to -4.0 in increments of 
0.1, corresponding to / — ► 1.26/ in the linear scale. 100 
simulated PSDs were used to determine the average model 
PSD. 10 4 randomly selected sets of PSDs were used to 
probe the Xdist distribution. 

Figure 7 shows the best fitting singly-broken PSD model 
shape and simulated data for each target, plotted in 
/ x Pf space. Figure 8 shows contour plots for the best- 
fitting singly-broken model fixed at the best-fitting value 
of 7. The presence of a constant Poisson power level dom- 
inates the steepest PSD slopes and leads to degeneracy in 
that most values of (3 above ~3 can lead to the same re- 
jection probability (e.g., upper limits on (3 for NGC 4151 
cannot be constrained); a relatively more minor effect is 
that excessive red-noise leak from very steep PSD slopes 
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(/? > 2.5) increases the errors, with the result that exces- 
sively large errors decrease the reliability of the values of 
the rejection probability. For each target, the best-fitting 
values of /?, 7, / c , and A and the likelihood of acceptance 
are summarized in Table 5. As before, errors on a single 
parameter arc la above the likelihood of acceptance for 
the best-fit value. The reader can also refer to the contour 
plots of Figure 8 for estimates of the absolute rejection 
probabilities for a given set of parameters. For Akn 564, 
no lower limit to 7 is given; models with blue-noise (7 < 0) 
PSD components are dominated by the aliased power, re- 
sulting in degenerate values of the rejection probability 
for any 7 < 0. The errors on A are determined from the 
RMS spread of the 10 4 randomly selected sets of simulated 
PSDs. 

For all six sources, the likelihood of acceptance improves 
when the break is added to the fit, as illustrated by the 
values of Act, the increase in likelihood of acceptance be- 
tween the unbroken and singly-broken power law fits, listed 
in col. (7) of Table 5. This in itself indicates problems 
with the pure power-law model. In four PSDs (Akn 564, 
NGC 3783, NGC 3516 and NGC 4151), the unbroken 
power law fit is accepted at less than 10% confidence while 
the singly-broken power law is accepted at greater than 
10% confidence. (That is, the rejection probability drops 
from greater than to less than 90%). These are signifi- 
cant improvements, corresponding to an increase in likeli- 
hood of acceptance of 0.9-2.9ct for these four sources. In 
the remaining two cases (Fairall 9 and NGC 5548) unbro- 
ken power law fits that were already reasonably acceptable 
(11.3% and 37.2% likelihood of acceptance) so the addi- 
tion of a break only improves the likelihood of acceptance 
by 0.4-0.6cr. 

4.3. A doubly-broken power law in NGC 3783? 

Inspection of Figure 7 suggests that the PSD NGC 3783 
is steep at high frequencies, then flattens out in / x Pf 
space (a slope of w — 1 in Pf space), then turns downward 
in / x Pf space (to a slope of « in Pf space) at the 
lowest frequencies probed. This behavior is reminiscent of 
low-state XRBs (e.g., Nowak et al. 1999a) and suggests an 
improvement upon the singly-broken fit in NGC 3783. 

In order to test for a second break, we employ a model 
of the form 

[ A h f < fi 

P(f)={ A h {f/f h )-\ f,<f<f h 
{ A h {f/f h )-0, f>f h 

where Ai is the PSD amplitude below the low-frequency 
break /;, where the PSD has zero slope. The intermediate 
slope, between the two break frequencies, is fixed at -1 in 
order to fix this model with the same number of free pa- 
rameters as the singly-broken model. Ah = A; x (fh/fi) -1 
and equals the PSD amplitude at the high frequency break 
fh] the PSD slope above fh is -[3. Break frequencies were 
tested in the log from -8.0 to -4.0 in increments of 0.1, cor- 
responding to / — » 1.26/ in the linear scale. The range 
of (3 tested was 1.1-2.4 in increments of 0.1. 100 simulated 
PSDs were used to determine the average model PSD. 10 4 
randomly selected sets of PSDs were used to probe the 
Xdist distribution. 

The best fitting model, plotted in Figure 9, has a low- 
frequency break at // = 2.00±f ^ x 10" 7 Hz, a high- 



frequency break at fh = 3. 98^47 x 10~ 6 Hz, and a power 
law slope of [3 = 2.0 ± 0.3 in Pf space above the high 
frequency break. The best-fitting power amplitudes are 
A h = 2.8t 1 ' 3 , x 10 3 Hz" 1 and Ai = 5.6±°j x 10 4 Hz" 1 . 
The likelihood of acceptance is 25.5% (rejection probabil- 
ity of 74.5%), only a modest improvement over the singly- 
broken power law model fit. The evidence in favor of two 
breaks in the NGC 3783 PSD is therefore only tentative 
at this point. 

One can speculate, however, that the value of 
7 = OTOig^g in the best-fitting singly-broken power law 
model for NGC 3783 actually reflects an average of the 
intermediate and white noise portions of a doubly-broken 
PSD. More significant, when 7 in the singly-broken model 
is fixed to the intermediate power law slope of -1 (similar 
to the best-fit singly-broken models for the other broad- 
line Scyfcrt 1 targets), the best fit, with [3 = 1-95^20 
and f c = 6.3ll3'gQxlO- 6 Hz, has a very low likelihood 
of acceptance, 2.2% (rejection probability of 97.8%). 

If Seyfert 1 PSDs do resemble XRB PSDs and have 
an an intermediate power law slope of -1, then the PSD 
of NGC 3783 does require a second break and resembles 
strongly the PSDs of low-state XRBs. 

4.4. Testing the stationarity of the light curves 

In order to combine the different time scale PSDs to 
synthesize a single broadband PSD, we must assume that 
the specific short time scale PSD measured in a single in- 
tensive XMM-Newton or Chandra observation is represen- 
tative of the average short time scale PSD that would have 
been measured if the source was monitored continuously 
for the entire ~ 3 yr period. We have tested this assump- 
tion using the five Chandra scans of NGC 3783, all taken 
with identical sampling intervals and durations. Each was 
used to measure separate short time scale PSDs, which 
have identical spectral coverage and suffer the same levels 
of red-noise leak and aliasing. This allows a straightfor- 
ward comparison, as shown in Figure 10. Note that all 
the PSDs show identical slopes to within the errors, and 
the amplitude normalization varies by ~40 percent. This 
small range of amplitudes is fully consistent with station- 
ary behavior, providing support for a key assumption of 
the PSD synthesis technique, and consistent with the lin- 
ear RMS-flux relation in the light curves of XRBs and 
Seyfert Is discussed in Uttley & McHardy (2001). 

The assumption of non-stationarity in the PSD applies 
not only to the short-term data, but must be applicable 
to the long-term light curves as well. To examine the as- 
sumption of stationarity on long time scales, we halved all 
the long-term light curves and applied the Monte Carlo 
analysis to the resulting PSDs by assuming the best-fitting 
model for each target. Reasonable fits were obtained for 
all targets, supporting stationarity of the PSDs on time 
scales of >1.5 years. 

5. DISCUSSION 

The fitting results are summarized in Figure 11, which 
allows the reader to directly assess the significance of fea- 
tures in the PSDs. These data are plotted in "model 
space" for ease of interpretation, with the data shifted rel- 
ative to the model as required to account for aliasing and 
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other distorting effects. This is the opposite of the earlier 
figures, done in "data space" in which the model is shifted 
to account for distortion effects. 

The plot shows clearly the inadequacy of describing 
these PSDs as traditional red-noise pure power-laws. All of 
the PSDs appear to show a similar departure from power- 
law behavior, a flattening to lower temporal frequencies, 
although the break strength and frequencies differ from 
object to object. The downturn is greatest in the PSDs of 
NGC 3783 (for which a "double break" is tentatively fa- 
vored) and Akn 564. NGC 3516 and NGC 4151 also show 
strong breaks. This is confirmed by statistical tests, which 
find that addition of a break improves the likelihood of ac- 
ceptance from less than 10% confidence to greater than 
10% confidence for all four of these PSDs (see § 4.2.). On 
the other hand, the PSDs of NGC 5548 and Fairall 9 are 
much steeper, with a turnover at lower temporal frequen- 
cies, and the addition of a break is not statistically sig- 
nificant. Only the upper limits to the break frequencies 
measured for Fairall 9 and NGC 5548 will be considered 
henceforth. We conclude that there is clear visual and sta- 
tistical evidence of a significant flattening in four of these 
objects, and visual indications of a break in the remaining 
two. The implications of these results are discussed below. 

5.1. Comparison with previous results 

The present analysis finds a break in the PSD of 
NGC 3516 at a frequency of 2 x 10~ 6 Hz (corresponding 
to a time scale of 6 days), a factor of ^5 higher frequency 
than the break found by EN99 in this object. The most 
likely cause for this discrepancy is the fact the EN99 forced 
the power law slope below the break to zero; the present 
analysis finds a much steeper low-frequency slope, lead- 
ing to the PSD break being modelled to occur at higher 
frequencies. 

In their analysis of the Akn 564 PSD, Pounds et al. 
(2001) also forced the low-frequency slope to zero and also 
found a power law slope above the break consistent with 
— 1, consistent with the model fit in this work. However, 
the present work finds a break frequency that is a factor of 
~2 higher; this is most likely due to the fact that Pounds et 
al. (2001) did not account for the large quantity of aliased 
power present in the PSD and allowed the relative normal- 
izations of individual PSD segments to vary. 

The sample of four objects in UMP02 was chosen on 
the basis of known strong X-ray variability, whereas the 
current sample was chosen for 2-10 keV flux (Piccinotti 
et al. 1982). The UMP02 sample could hence be biased 
towards rapid variability. Two objects in the UMP02 sam- 
ple, NGC 3516 and NGC 5548, overlap with ours, and the 
PSD analyses yield consistent results. The two other ob- 
jects, NGC 5506 and MCG-6-30-15, are modelled to show 
very high frequency breaks (^5xlO _5 Hz) when a break to 
a fixed low-frequency slope of -1 is assumed; however it 
also should be noted that fitting these objects' broadband 
PSDs with a model incorporating a fixed low-frequency 
slope of yields PSD breaks located an order of magnitude 
lower in frequency. NGC 5506, from the UMP02 sample, 
is a Seyfert 2 and will not be discussed further. MCG-6- 
30-15 will be included along with the other six Seyfert 1 
galaxies in the ensuing subsections below. 

In their PSD analysis, Hayashida et al. (1998) defined 



a characteristic variability frequency in the PSDs of AGN 
and Cyg X-l as being the frequency at which the / x Pf 
power crossed a certain threshold value (10~ 3 ). For the 
PSD shapes observed in XRBs and Seyferts, such a fre- 
quency is dependent on the PSD break frequency, high- 
frequency slope and normalisation. However, all these 
quantities can vary in XRBs, leading to degenerate val- 
ues of the characteristic variability frequency for a given 
black hole mass. It is also more plausible that the break 
frequencies, and not the high-frequency slopes, are more 
closely connected to the black hole mass (see § 5.3). Fur- 
thermore, the distortion effects of red-noise leak and alias- 
ing were not taken into account when determining those 
PSDs. 

Complementary to the work of Pounds et al. (2001), 
Papadakis et al. (2002) claim evidence for a break in PSD 
slope from <~ —1.3 to ~ —1.7 above a temporal frequency 
of <~ 2 x 10~ 3 Hz in Akn 564. However, Papadakis et 
al. (2002) estimate the high-frequency noise level simply 
from the average flux errors. Consequently, any small er- 
ror in this noise level determination, which dominates the 
highest frequencies probed, requires a subtantial change in 
high-frequency PSD slope. Additionally, that ASCA data 
only were used leaves about 2/3 of a order of magnitude 
of temporal frequency uncovered, due to Earth-occultation 
gaps; more importantly, distortion effects are not consid- 
ered. 

5.2. Comparison with X-ray binaries 

The Seyfert PSDs are much more poorly defined than 
those seen in low-state XRBs (e.g., Nowak et al. 1999a, 
Nowak et al. 1999b). McHardy (1988) first made the link 
between the two types of compact systems based on similar 
slopes of the EXOSAT power-law PSDs (see also Lawrence 
& Papadakis 1993). The detection of Seyfert 1 PSD breaks 
(eg, EN99), on time scales that appear to scale linearly 
with the breaks seen in XRBs, helped confirm this picture. 
Furthermore, the energy spectra of both types of objects 
are similar, featuring a non-thermal hard X-ray coronal 
power law which steepens as the source brightens (e.g., 
Markowitz & Edelson 2001, Wilms et al. 1999), a 6.4 keV 
Fe Ka fluorescence line modelled to orignate in the regime 
of strong gravity (Tanaka et al. 1995, Fabian et al. 1995, 
Wilms ct al. 1999), and a Compton reflection hump around 
20-30 keV. These spectral components, along with the 
rapid X-ray variability and PSD similarities such as and 
the energy dependence of the PSDs (Nandra & Papadakis 
2001) and apparent phase lags (Papadakis, Nandra & 
Kazanas 2001), supporting common X-ray emission and 
variability mechanisms in both classes of compact objects. 

The observed PSD break frequencies in Seyferts are 
about 6-7 orders of magnitude smaller than those see in 
low-state XRB PSDs. It is remarkable to note that this 
ratio is approximately the same as that of the X-ray lumi- 
nosities and putative black hole masses (10 Mq for Cyg X- 
1, Herrero et al. 1995; 1O 7 ~ 8 M for Seyferts, e.g., Kaspi 
et al. 2000). This scaling supports the picture of similar 
variability processes operating in both types of accreting 
compact systems. It is noted that this type of scaling is 
subject to the caveat of Seyferts and low-state XRBs differ- 
ing only in mass, physical size, luminosity and variability 
time scale. In fact, differences between the two classes of 



8 



objects, most notably lower inner accretion disk tempera- 
tures in Scyferts than seen in low-state XRBs, prevent the 
analogy from being solidified. 

However, key questions remain in this picture. Are the 
PSD shapes really identical between the two types of sys- 
tems? The fact that XRB PSDs do not strictly adhere 
to a doubly-broken power law (e.g., GX 339-4, Nowak 
et al. 1999b) and that our analysis only quantifies the 
overall broad PSD shape in terms of simple broken power 
laws hinder this issue. The sharpness of these models' 
breaks lacks a physical basis, but the data are not ad- 
equate enough to accurately constrain the 'bluntness' of 
the break. Do the Seyfert PSD breaks seen correspond 
to the high or low frequency breaks seen in XRB PSDs? 
The tentative finding of two PSD breaks in NGC 3783 
matches closely the double breaks in XRB PSDs, and the 
ratio of the break frequencies in both objects are nearly 
the same (~ 20 for NGC 3783, -20 for Cyg X-l by Nowak 
et al. 1999a; 25 for the ratio of the average Cyg X-l break 
frequencies, Belloni & Hasinger 1990), though the low- 
frequency break in Cyg X-l is somewhat more variable 
and so these ratios should be taken with a grain of salt. 
For the other Seyferts, however, this issue is not as clear. 
The singly-broken model fits for the PSDs of NGC 3516, 
NGC 4151, and to a lesser extent, NGC 5548 and Fairall 9, 
appear similar to the high frequency break in XRBs in 
terms of incorporating a break above which the power law 
slope is —-2 and below which the slope is ~-l (see UMP02 
for additional discussion). In contrast to the broad-line 
targets, the best-fitting model for the soft-spectrum source 
Akn 564 seems to mimic the low-frequency breaks seen in 
low-state XRBs, that from the intermediate slope of -1 
above the break to zero below it (consistent in shape to 
the low-frequency break in the doubly-broken model fit to 
the NGC 3783 PSD). Alternatively, we caution the reader 
that soft-spectrum targets like Akn 564 may in fact ex- 
hibit fundamentally different variability characteristics in 
the 2-10 keV band from broad-line Seyfert 1 galaxies (see 
e.g., Edelson et al. 2002), and the entire underlying broad- 
band PSD shape may in fact be significantly different in 
shape from low-state XRBs. The limited dynamic range of 
the present PSD for this target hinders this question. The 
low-frequency PSD flattening to zero slope in Akn 564 
presents difficulty in solidifying the inviting analogy be- 
tween soft-spectrum Seyfert 1 galaxies such as Akn 564 
and high-state XRBs (which show 1// noise down to fre- 
quencies of a few x 10~ 3 Hz, e.g., Churazov, Gilfanov & 
Revnivstev 2001), as both types of systems are suspected 
of accreting at relatively high fractions of the Eddington 
limit (e.g., Pounds, Done & Osborne 1995) and both show 
a steep energy spectrum (r — 2.4, e.g., Vaughan et al. 
1999; Lcighly 1999). Finally, what physical process is re- 
sponsible for the variability? There are few specific models 
available, though some are discussed briefly in § 5.4. 

5.3. Correlations between break time scale, luminosity 
and mass estimates 

For the seven Seyfert Is under consideration, Table 6 
lists bolometric luminosity Lboi, black hole mass esti- 
mate Mbh, and PSD break time scale Todays) from the 
singly-broken power law model fits (for NGC 3783, the 
high-frequency break from the doubly-broken model fit is 



used). The time scale for MCG-6-30-15 is from the 'high- 
frequency break' model used by UMP02. Bolometric lu- 
minosity is calculated using the mean bolometric correc- 
tion of Padovani & Rafanelli (1988), Lboi = 27i 2 -io koV- 
It should be noted that for the soft-spectrum Seyfert 1 
Akn 564, the 2-10 keV X-ray to bolometric luminosity con- 
version of Padovani & Rafanelli (1988) should be treated 
as a conservative lower limit due to the large EUV ex- 
cesses present in soft-spectrum Seyfert Is (see, e.g., Turner 
et al. 2002). The bolometric luminosity for MCG-6-30-15 
is taken from Reynolds et al. (1997). All reverberation- 
mapped masses are taken from Kaspi et al. (2000), except 
for NGC 3516, which is from Wanders & Home (1994). 
Akn 564 does not have a highly reliable reverberation- 
mapped mass estimate; the upper virial mass estimate is 
< 8 x 1O 6 M by Collier et al. (2001), based on the lag 
between UV continuum flux and Lya A1316 variations. 
Reverberation mapping has not yet been done for MCG- 
6-30-15, but its black hole mass estimate may be derived 
as follows: Reynolds (2000) noted that the bulge lumi- 
nosity of the host yields an approximate bulge mass of 
3 x 10 9 M Q . The FWHM H/3 is 1700 km/sec (Pineda et 
al. 1980), which, along with its steep energy index nearly 
makes it consistent with the properties of most 'officially'- 
classified NLSyls such as Akn 564. Using the improved 
black hole mass- bulge relations from Wandel (2002; see 
their Figures 3 & 4), for which NLSyl black hole masses 
lie about an order of magnitude below those for broad-line 
AGN and normal galaxies, the black hole mass for MCG- 
6-30-15 is likely to be of order - 1 x 10 6 M Q , the value 
adopted henceforth in this paper. 

Figure 12 shows Lboi and black hole mass plotted versus 
time scale. For luminosity versus time scale (top panel), 
the Pearson correlation coefficient r is 0.49, with a prob- 
ability P r of obtaining that value of r by chance of 0.27 
for seven points. The bottom panel of Figure 12 shows 
that there is an apparent correlation between mass and 
time scale, with r = 0.79, and P r — 3.3 x 10~ 2 for seven 
points. It should be noted that for NGC 3783, the low fre- 
quency PSD break from the doubly-broken model fit is not 
consistent with the above correlation. When the singly- 
broken model time scale is used for NGC 3783 (as op- 
posed to the high-frequency break from the doubly-broken 
fit), r decreases to 0.48 (P r — 0.27 for seven points) in 
the luminosity-time scale plot, and r decreases to 0.78 
(P r = 3.9 x 10~ 2 for seven points) in the mass-time scale 
plot. 

A more significant difference arises when plotting 
£2-10 koV against time scale, using the X-ray luminos- 
ity for MCG-6-30-15 from UMP02, where r drops to 0.27 
(P r = 0.56 for seven points). Finally, because the PSD 
break for MCG-6-30-15 was obtained from a fit by UMP02 
where the low frequency slope was fixed to -1 and fixing 
the low frequency to results in a break time scale of 2.26 
days, we note that using this time scale instead results in 
r = 0.52 (P r = 0.23 for seven points) for the luminosity- 
time scale plot and r = 0.62 (P r = 0.14 for seven points) 
for the mass-time scale plot. 

The apparent correlation between mass and time scale 
seen in Figure 12, stronger than the luminosity-time scale 
relation, supports the validity of the reverberation map- 
ping method of black hole mass estimation; furthermore, 
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one could speculate that mass may be more relevant than 
luminosity in governing a given Seyfert's X-ray timing 
properties. The mass-time scale correlation is consistent 
with the linear relation T( days ) = Mbh/10 6 ' 5 M q , but this 
is certainly not conclusive with only seven objects span- 
ning a narrow range of luminosity and black hole mass. 
The mass-time scale correlation and the similar shapes of 
the PSD breaks are consistent with a picture in which the 
PSDs of all the Seyfert Is considered here have a similar, 
universal broadband PSD shape (that of low-state XRBs) 
in which the PSD breaks move towards lower temporal 
frequency as black hole mass increases. This is consistent 
with the expectation that compact accreting black hole 
systems with a relative larger black hole mass and larger 
Schwartzschild radius, i?sch, have a a larger X-ray emis- 
sion region, requiring a comparatively longer duration for 
the source to achieve a given amplitude of flux variabil- 
ity and longer characterstic variability time scales. This 
paradigm is manifested in findings of an anticorrclation be- 
tween X-ray luminosity and excess variance, as measured 
over a fixed duration (Nandra et al. 1997, Turner et al. 
1999, Markowitz & Edelson 2001). 

It is also interesting to note that if one extrapolates the 
linear mass-time scale relation down to Cyg X-l's mass of 
10 M Q , the inverse of the resulting time scale is 3.7 Hz, 
remarkably close to Cyg X-l's mean high-frequency break 
of 3.3 Hz, again supporting the notion of a common vari- 
ability mechanism in both types of objects. This is also 
suggestive that there is only a small departure from linear- 
ity of the mass-time scale relation over 6-7 decades. Fur- 
thermore, this suggests that the PSD breaks detected are 
high-frequency breaks, supported by the similar changes in 
slope from ~-2 to ~-l (with Akn 564 being an exception). 

Finally, we note that the claim of a high-frequency break 
at - 2 x 10~ 3 Hz in the PSD of Akn 564 by Papadakis 
et al. (2002) does not fall anywhere near this mass-time 
scale relation; Papadakis et al. (2002) note this fact, and 
consider the possibility that the break detected may not 
be indicative of black hole mass. 

5.4. Correlations between PSD amplitude, luminosity, 
and mass estimates 

Figure 13 shows the bolometric luminosity and black 
hole mass plotted against A, the PSD amplitude at the 
break (in Pf space) in the singly-broken model fits (high- 
frequency break from the doubly-broken fit for NGC 3783). 
Also included is MCG-6-30-15's break amplitude from 
UMP02, A = 213 ± 8 Hz" 1 . Under the assumption 
that it is a 'low-frequency' break, that of Akn 564 is 
excluded. Though there are are only six data points, 
there appears to be strong correlation between black hole 
mass and PSD break amplitude, with Pearson r = 0.95 
(P r = 3.1 x 10~ 3 ). There is a moderate correlation be- 
tween iboi and PSD break amplitude, with r = 0.72 
{P r = 0.10). 

These correlations reflect how the PSD, still assumed to 
have the same universal shape as low-state XRBs, not only 
moves downward in temporal frequency but also upward 
in power at the break frequency as black hole mass in- 
creases. In a sample of nine broad-line Scyferts (including 
all six broad-line Seyferts discussed thus far), Markowitz 
& Edelson (2001) saw that the dependence of long-term 



excess variance (measured over 300 d) on luminosity was 
much shallower than the dependence on short-term (Id) 
excess variance. This indicated that the long-term excess 
variances likely probed down to temporal frequencies be- 
low targets' respective PSD breaks. As a given PSD moves 
towards lower temporal frequencies as black hole mass in- 
creases, the power level at the break must increase so that 
long-term excess variance stays similar. The above corre- 
lation arises mainly because the peak power in / x Pf 
space is broadly the same for all targets (/ x Pf ~ 0.01), 
similar to that of Cyg X-l, as discussd by e.g., UMP02, not 
only strengthening the notion of similar variability mech- 
anisms in Seyferts and XRBs, but also implying similar 
numbers of varying X-ray emitting regions. 

5.5. Testing physical time scales 

Since the observed cutoff is expected to relate to the 
physical process that generates the variability, one can ex- 
plore the relevance of various physical processes to the 
origin of the X-ray variability by comparing physical time 
scale predictions to the measured PSD time scales. EN99 
discuss the relevance of light-travel time effects and the 
orbital time scale for a standard a-disk (Shakura & Sun- 
yaev 1973), comparing them to the PSD break observed in 
NGC 3516. EN99 note that for NGC 3516, and we note for 
our Seyferts here, that these time scales are both generally 
too short to be directly associated with the days-to-weeks 
PSD time scales. Additionally, the radial drift /viscous 
time scale is usually a few orders of magnitude longer 
than observed PSD time scales. However, the thermal and 
acoustic variation time scales for a thin disk (days to weeks 
given the putative black hole masses considered; Maraschi, 
Molendi & Stella 1992; Treves, Maraschi & Abramowicz 
1988) are closest to the measured characterstic variability 
time scales. We hence speculate that the mechanism of 
X-ray variability may be tied to physical processes which 
incorporate thermal or acoustic thin-disk variations. If 
this is indeed the case, then the characteristic time scales 
and black hole mass estimates can be used to constrain the 
X-ray emitting location. Assuming reasonable values for 
the scale height and viscosity, and assuming that the thin 
disk approximation holds over the entire radius range, all 
of the targets' acoustic and thermal time scales are con- 
sistent with the variable X-ray emission originating within 

^30i?Sch- 

However, more specific models can be considered. For 
instance, in the pulse avalanche model of Poutanen & 
Fabian (1999), the variable X-ray emission is produced in 
an inhomogeneous, stochastic system of magnetic flares 
inflating and detaching above the inner accretion disk, 
producing X-ray flares via inverse-Compton scattering of 
softer photons, but with events occurring in correlated 
avalanches. The resulting light curves can yield doubly 
broken PSDs similar to those of low-state XRBs for ap- 
propriately tuned parameters, with a break between 
noise and f" 1 noise being associated with the longest time 
scale (days to weeks for the broad-line Seyfert Is consid- 
ered here) of an individual flare. A low- frequency PSD 
break from / _1 to white noise is associated with the dura- 
tion of the avalanches. However, predicted flare durations 
are generally too short to be directly associated with PSD 
breaks. 
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Along another vein, the X-ray variability in the model of 
Lyubarskii (1997), is caused by variations in the accretion 
flow (e.g., due to turbulence) propagating from large radii 
inwards to the location of the X-ray corona. Such fluc- 
tuations in viscosity are dependent on the local viscous 
time scale, which decreases with decreasing radius. For 
standard thin disks, the viscous timescale is too long to be 
associated with the measured PSD breaks. However, if the 
accretion flow is geometrically thick, such as an advection- 
dominated accretion flow (ADAF; Narayan & Yi 1994) or a 
radiation-supported thick disk (e.g. see Abramowicz 1988, 
Treves, Maraschi & Abramowicz 1988) the viscous time- 
scale tends toward the thermal time-scale. Therefore the 
PSD breaks seen might plausibly correspond to a viscous 
time-scale in an ADAF or radiation-supported thick disk. 

6. CONCLUSIONS 

We have systematically constructed high-quality broad- 
band PSDs for a sample of six Seyfert 1 galaxies. These 
high-quality PSDs cover exceptional dynamic ranges, con- 
tinuously spanning up to or beyond 3-4 orders of mag- 
nitude in temporal frequency. We use the Monte Carlo 
technique of UMP02 to determine adequate errors on each 
binned PSD point and to account for distortion effects, 
and we characterize the underlying PSD shape to look for 
low-frequency flattening below a characteristic break fre- 
quency. Four targets (Akn 564, NGC 3783, NGC 3516 and 
NGC 4151) show significant evidence for low-frequency 
flattening, with characteristic variability time scales in the 
range of several to tens of days. For the two most mas- 
sive and luminous targets in the sample (Fairall 9 and 
NGC 5548), we expect that continued long-term monitor- 
ing and extending the PSD to probe lower frequencies will 
confirm their break frequencies, which are the lowest in 
the sample. The low-frequency flattening seen in these 
objects is remarkably similar to what is seen in low-state 
X-ray binary PSDs, strengthening the argument that sim- 
ilar emission processes occur in both types of compact ac- 
creting systems, spanning a factor of ^10 6-7 in luminosity 
and putative black hole mass. 

All of the PSDs studied are consistent in shape with at 
least portions of low-state XRB PSDs. The finding of two 
breaks in the PSD of NGC 3783 is only tentative, but it 
still heralds a new era in Seyfert PSD analysis: those PSDs 
derived from EXOSAT data generally lacked the tempo- 
ral frequency range to reveal even a single break, and the 
access to long time scales afforded by RXTE allowed detec- 
tion of a single break (EN99, Pounds et al. 2001, UMP02, 
this work). However, as double power law breaks are rou- 
tinely detected in XRB PSDs, only Seyfert PSDs which 
also show two breaks will exhibit the most direct and con- 
vincing evidence for a link between the two classes of sys- 
tems. More long-term monitoring is therefore crucial for 
locating low-frequency PSD breaks. 

The PSD break frequencies detected are model- 
dependent, but when compared to reverberation-mapped 
black hole masses, they are consistent with a linear mass- 
time scale relation that extends not only throughout the 
Seyfert range, but down to the mass and PSD time scale 
of Cyg X-l, providing even further support for similar 
variability processes in Seyferts and XRBs. The corre- 

11 See http://www.src.le.ac.uk/lobstcr 



lation, however, is hampered by small-number statistics, 
and more accurate black hole mass estimates for Akn 564 
and MCG-6-30-15 are needed to clarify the situation; this 
is especially important for MCG-6-30-15 since its mass es- 
timate is very crude, yet it is used to anchor the low-mass 
end of the black hole mass-time scale relation shown in 
Figure 12. 

It is encouraging to note that this agreement between 
X-ray variability time scale and reverberation mapping 
arises despite: numerous caveats and assumptions in both 
analysis methods (e.g., the assumption of broad-line emis- 
sion clouds being on Keplerian orbits, Peterson & Wandel 
1999), each of which is expected to probe different regions 
of the central engine. 

Finally, we draw the reader's attention back to Fig- 
ure 11, as we believe such diagrams will play a large role 
in future Seyfert 1 PSD analysis. It is a plot of / x Pf as 
a function of / in model space, while most previous PSD 
plots were of Pf as a function of / in data space. In this re- 
gard, it mimics the "unfolded" (model-space) v x F v plots 
of early X-ray energy spectral analyses. That was done 
to give the clearest possible visual picture of the instrinsic 
energy spectrum itself, undistorted by the vagaries of the 
response matrix of a particular detector. Likewise, this 
plot gives the clearest possible picture of the underlying 
PSD, after removing the distortions unique to the partic- 
ular sampling pattern employed. Further, just as a peak in 
the v x F v plots indicated a region of maximum luminosity 
per decade of frequency, so a well-defined peak in the / x Pf 
PSD plot occurs in the temporal frequency band at which 
most of the variability power is produced. These data in- 
dicate that the "characteristic variability time scales" of 
these Seyfert Is are typically of order a few days. 

Establishing that breaks are a systematic feature of their 
PSDs represents a milestone in the description of Seyfert 1 
X-ray variability. Here too there is analogy with the early 
energy spectral studies: the low-resolution energy spec- 
tra from early missions (e.g., Ariel V) were wcll-modclcd 
as "canonical" T ~ 1.7 power-laws, but a critical break- 
through came when Wilkes & Elvis (1987) unambiguously 
established the existence of a significant departure, the 
"soft excess," from the Einstein data. As the spectral res- 
olution and bands improved over the following 15 years, 
many more spectral features have been found (e.g., narrow 
absorption lines that make up the "warm absorber," the 
iron K Q fluorescence line, the Compton reflection compo- 
nent). Now this work has clearly established the perva- 
siveness of breaks in Seyfert 1 PSDs that were modeled as 
simple power-laws in previous, lower dynamic range PSDs. 
We can not be certain that this is strictly a broken power- 
law; many other paramcterizations are of course possible. 
We can be fairly confident, however, that the current de- 
scription of Seyfert 1 PSDs will be refined and eventually 
overturned, just as the "canonical" energy spectral power- 
law was 15 years ago, as we further probe long time scales 
with RXTE and short time scales with XMM-Newton, and 
the flight of Lobster-ISS 11 starting in 2009 increases the 
number of Seyfert 1 PSDs more than tenfold. 
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Table 1 

Source Parameters for the Seyfert Galaxies 



Source 
Name 


^2-lOfceV 

(loglO(erg s- 1 )) 


z 


FWHM H/3 
(km sec^ 1 ) 


I^-IO 


Fairall 9 


43.97 


0.047 


5900 


2.18 


NGC 5548 


43.50 


0.017 


5610 


1.79 


Akn 564 


43.49 


0.024 


720 


2.68 


NGC 3783 


43.22 


0.010 


2980 


1.77 


NGC 3516 


42.86 


0.009 


6800 


1.60 


NGC 4151 


42.62 


0.003 


5000 


1.65 



Note. — Targets are ranked by unabsorbed 2-10 keV luminosity, listed in 
Col. (2). Luminosities were calculated with the online W3PIMMS tool and 
using the long-term mean RXTE count rates, X-ray photon indices, assum- 
ing a cold absorption column equal to the Galactic column (for NGC 4151, a 
column density 230 x the Galactic column was used; see Schurch & Warwick 
2002), and calculated assuming H = 75 km s _1 Mpc -1 and q = 0.5. Red- 
shifts (col. [3]) were obtained from the NED database. H/3 values (col. [4]) 
were taken from Turner ct al. (1999); see references therein. The photon 
indicies in col. (5) were obtained from the Tartarus database, derived from 
a simple power law fit covering the 2-10 keV bandpass, except for Akn 564, 
from a spectral fit by Vaughan et al. (1999), NGC 3783, from a spectral fit 
to Chandra data by Kaspi et al. (2001), and NGC 4151, from a spectral fit 
to XMM-Newton data by Schurch & Warwick (2002). 



Table 2 
Sampling Parameters 



Source 


Time 




MJD Date 






Fraction 






Name 


Scale 


Instrument 


Range 


AT 


Npts 


missing (%) 




(%) 


Fairall 9 


Long 


RXTE 


51180.59-52311.18 


4.27 d 


266 


7.9 


1.54 


39.2 ± 1.8 


Fairall 9 


Medium 


RXTE 


52144.88-52178.97 


3.2 hr 


258 


13.2 


2.09 


14.7 ± 0.7 


NGC 5548 


Long 


RXTE 


50208.07-52417.02 


14 d 


158 


0.0 


4.82 


34.5 ± 0.9 


NGC 5548 


Medium 


RXTE 


52091.70-52125.40 


3.2 hr 


255 


8.2 


4.51 


26.3 ± 1.2 


NGC 5548 


Short 


XMM 


52097.68-52098.63 


900 s 


93 


0.0 


2.99 


5.3 ± 0.4 


Akn 564 


Long 


RXTE 


51179.58-52310.19 


4.27 d 


266 


7.9 


1.73 


29.2 ± 1.3 


Akn 564 


Medium 


RXTE 


51694.86-51726.51 


3.2 hr 


239 


13.0 


2.00 


34.9 ± 1.7 


NGC 3783 


Long 


RXTE 


51180.55-52311.25 


4.27 d 


266 


9.4 


7.09 


22.3 ± 1.0 


NGC 3783 


Medium 


RXTE 


51960.17-51980.05 


3.2 hr 


151 


7.9 


5.65 


12.8 ± 0.8 


NGC 3783 


Short 


Chandra 


51964.80-51966.74 


2000 s 


84 


0.0 


0.39 


12.3 ± 0.9 


NGC 3783 


Short 


Chandra 


51967.40-51969.38 


2000 s 


84 


0.0 


0.39 


11.9 ± 0.9 


NGC 3783 


Short 


Chandra 


51978.04-51979.94 


2000 s 


84 


0.0 


0.41 


11.5 ± 0.9 


NGC 3783 


Short 


Chandra 


51999.17-52001.11 


2000 s 


84 


0.0 


0.58 


15.1 ± 1.2 


NGC 3783 


Short 


Chandra 


52086.43-52088.38 


2000 s 


84 


0.0 


0.49 


13.5 ± 1.0 


NGC 3516 


Long 


RXTE 


50523.03-51593.40 


4.27 d 


252 


15.1 


3.77 


37.8 ± 1.8 


NGC 3516 


Medium 


RXTE 


50523.03-50659.09 


12.8 hr 


256 


22.7 


4.21 


28.0 ± 1.4 


NGC 3516 


Short 


RXTE 


50590.01-50594.22 


1200 s 


304 


15.1 


3.85 


7.6 ± 0.3 


NGC 3516 


Short 


RXTE 


50916.34-50919.67 


1200 s 


241 


29.1 


5.09 


10.4 ± 0.6 


NGC 4151 


Long 


RXTE 


51179.56-51964.65 


4.27 d 


185 


6.0 


14.86 


38.1 ± 2.0 


NGC 4151 


Medium 


RXTE 


51870.60-51904.79 


3.2 hr 


259 


11.2 


7.43 


18.4 ± 0.9 


NGC 4151 


Short 


XMM 


51900.48-51901.14 


1100 s 


53 


0.0 


3.04 


7.8 ± 0.8 



Note. — Targets are ranked by 2-10 keV luminosity. Col. (5), AT samp , is the sampling interval. Col. (8), fi, is the mean 
count rate. RXTE count rates are per 1 PCU. 
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Table 3 

PSD Measurement Parameters 



Source 


1 lme 


£ 

Jmin 


£ 

JNyq 


Temp. Freq. Range 


rpsn 


Name 


bcale 


(Hz) 


(Hz) 


Spanned (decades) 


(Hz ) 


Fairall 9 


Long 

o 


1.0 x 10~ 8 


1.4 x 10~ b 


3.6 


1360 


Fairall 9 


Medium 


3.4 x 10~ 7 


4.4 x 10" 5 




30.9 


NGC 5548 


LonE 


5.2 x 10" y 


1.4 x 10~ b 


5.0 


577 


NGC 5548 


Medium 


3.4 x 10~ 7 


4.4 x 10" 5 




9.4 


NGC 5548 


Short 


1.2 x 10" 5 


5.6 x 10" 4 




0.69 


Akn 564 


LonE 


1.0 x 10~ 8 


1.4 x 10~ b 


3.6 


1050 


Akn 564 


Medium 


3.7 x 10~ 7 


4.4 x 10~ 5 




33.8 


NGC 3783 


Long 


1.0 x 10" 8 


1.4 x 10" b 


4.4 


1210 


NGC 3783 


Medium 


5.8 x 10~ 7 


4.4 x 10~ 5 




6.87 


NGC 3783 


Short 


5.9 x 10~ 6 


2.5 x 10~ 4 




5.12 


NGC 3783 


Short 


5.9 x 10~ 6 


2.5 x 10~ 4 




5.14 


NGC 3783 


Short 


5.9 x 10~ 6 


2.5 x 10~ 4 




4.94 


NGC 3783 


Short 


5.9 x 10~ 6 


2.5 x 10~ 4 




3.45 


NGC 3783 


Short 


5.9 x 10" 6 


2.5 x 10" 4 




4.13 


NGC 3516 


Long 


1.1 x 10" 8 


1.4 x 10" b 


4.6 


177 


NGC 3516 


Medium 


8.5 x 10~ 8 


1.1 x 10" 5 




13.8 


NGC 3516 


Short 


2.8 x 10~ 6 


4.2 x 10~ 4 




0.37 


NGC 3516 


Short 


3.5 x 10" 6 


4.2 x 10" 4 




0.23 


NGC 4151 


Long 


1.1 x 10" 8 


1.4 x 10" b 


4.6 


44.2 


NGC 4151 


Medium 


3.4 x 10~ 7 


4.4 x 10" 5 




4.7 


NGC 4151 


Short 


1.2 x 10" 5 


4.5 x 10~ 4 




0.67 



Table 4 

Results For Unbroken Power Law Model Fit 



Best- Best-Fitting A a Likelihood 
Target Fitting f3 (units of Hz -1 ) of acceptance (%) 



Fairall 9 1.60to^ 2.7±0.4 x 10 3 11.3 

NGC 5548 1.65±°1o 3.4±g; \ x 10 3 37.2 

Akn 564 0.75to' 4 o 4.5±0.5 x 10 3 0.6 

NGC 3783 1.25 5.5±0.5 x 10 3 <0.01 

NGC 3516 1.80±g;|g 9-3±i;| x 10 3 6.6 

NGC 4151 1.90±g;|g 7.8^ : 4 x 10 3 1.8 



Note. — Results from fitting the PSDs with an unbroken power law model. 
A a is the power amplitude at / = 10~ 6 Hz. The likelihood of acceptance is 
defined as (1 - R), where R is the rejection confidence. No error is assigned for 
the best-fitting /? for NGC 3783. Upper limits on (3 for Fairall 9, NGC 3516 and 
NGC 4151 cannot be constrained due to a red-noise leak bias discussed in the 
text. 
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Table 5 

Results For Singly-Broken Power Law Model Fits 













Likelihood 


Target 


7 


P 


fc (Hz) 


A (Hz- 1 ) 


of acceptance (%) Act 



Fairall 9 


1 m +11 ° 

i - iu -0.60 


o on+ - 65 

z - zu -0.20 


3.98+2^ x 10" 


-7 


5.4t°- 8 7 x 10 4 


23.5 


0.6 


NGC 5548 


1 1 c+0.50 
1 - lo -0.65 


2 r)5+°- 80 


6.31^x10- 


-7 


2.5±g"5 x io 4 


87.4 


0.4 


Akn 564 




1 20+ - 25 


1.59^'^xlO- 


-6 


1.9^; 2 x 10 4 


97.3 


2.9 


NGC 3783 




i 9n +0 15 


2.ooig-^xio- 


-6 


1.0±0.1 x 10 4 


12.7 


2.0 


NGC 3516 


1 io+°- 40 

1 - lu -0.30 


2 oo+°- 55 


2.00^'^xlO- 


-6 


7.9±g;| x 10 3 


81.1 


1.9 


NGC 4151 


1 10+ 1 - 35 


2 in +1 - 90 


1.26t i: ^xl0- 


■6 


1.3±0.2 x 10 4 


10.8 


0.9 



Note. — Results from fitting the PSDs with the singly-broken power law model. 7 is the low frequency 
power law slope, (3 is the high-frequency power law slope, and A is the power amplitude at the break 
frequency f c . The likelihood of acceptance is again defined as (1 - R), where R is the rejection confidence. 
Act quantifies the increase in likelihood of acceptance between the unbroken and singly-broken power law 
fits. 



Table 6 

Time Scale, Bolometric Luminosity, and Black Hole Mass Estimate 





PSD Break 


Lbol 


B.H. Mass 


Target 


Time Scale (days) 


log(erg s" 1 ) 


(1O 7 M ) 


Fairall 9 


>18.3 


45.61 


8-3±t;| 


NGC 5548 


>4.6 


45.01 


9-4±i;I 


Akn 564 


7 3+ 11 

' - J -5.5 


44.81 


0.8 


NGC 3783 


2.9±U 


44.65 


1 10+ 107 


NGC 3516 


5.8+H 


44.36 


2.0±0.3 


MCG-6-30-15 


23+ - 67 


44.24 


0.1 


NGC 4151 


q 2+ 37 

y - Z -5.5 


44.06 


.70 



Note. - 



Targets are ranked by Lb i- See text for details. 
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Fig. 6. — Comparison of the best-fitting unbroken power law model against the observed PSDs, plotted in / X Pf space. The observed 
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Fig. 7. — Comparison of the best-fitting singly-broken power law model fits against the observed PSDs, plotted in / X Pf space. All lines 
and symbols are the same as in Figure 6. 
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Fig. 8. — Contour plots showing the errors on the best-fitting singly-broken PSD model shape. For each target, a slice through the 
three-dimensional fit parameter space at the best-fitting low-frequency slope is shown. Solid lines indicate the 99%, 95%, and 68% rejection 
probability levels. Lighter shading denotes a high rejection probability; darker shading denotes a low rejection probability. 
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Fig. 9. — Comparison of the best-fitting doubly-broken power law model fit to the observed NGC 3783 PSD, plotted in / X Pf space. All 
lines and symbols are the same as in Figure 6. 
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Fig. 10. — The five NGC 3783 Chandra short-term binned PSDs. From the first-observed to last observed, the PSDs arc marked by filled 
circles, crosses, open squares, filled stars, and filled triangles. The constant level of power due to Poisson noise has not been subtracted, though 
the Poisson noise levels are very similar for each observation. PSD errors are omitted here; they are determined later for a given PSD model 
shape. The best-fit lines (from first- to last-observed, respresented by bold solid, bold dotted, solid, dashed, and dotted lines respectively) 
were calculated assuming equal weighting to all points, but they are meant only to guide the eye, and are not an accurate representation of 
the intrinsic, underlying PSD slope due to the presence of red-noise leak and power due to Poisson noise. 



25 



tttt| — i i i iiiii| — i i i iiiii| — i i 1 1 1 1 ii| — i i 1 1 1 u i| — i i i rrrpnT] — i i i iiiii| — i i i iiiii| — i i i iiiii| — i i i iiiii| — i i rrrrpTTT] — i i i iiiii| — i i i iiiii| — i i i n iii| — i i i iiiii| — rr 



10 5 y 



ml 



NGC 3783 

ill I I 




■"lull 



NGC 3516 

ill I I 



iTTii ill 



NGC 4151 

ill I I 



ill 3 



10" 8 10~ 7 10" 6 10" 5 10~ 4 10~ 8 10~ 7 10~ 6 10" 5 10~ 4 10~ 8 10~ 7 10" 6 10~ 5 10~ 4 



Temporal Freq. (Hz) 



Temporal Freq. (Hz) 



Temporal Freq. (Hz) 



Fig. 11. — The best-fitting model fits (doubly-broken power law model lor NGC 3783, singly-broken power law model for all other targets) 
are shown with the effects of Poisson noise, aliasing and red noise leak subtracted off. The solid line indicates the underlying, intrinsic PSD 
model shape. Symbols represent the differences between the average distorted model and the observed PSD data, plotted relative to the 
underlying PSD model shape. 



26 



CO 



CD 



o 
PQ 



10 45 b 



10 44 b 







10 8 b- 



co 

CO 

« io 7 b- 

O 

o 

i — i 
PQ 

10 6 



> 



=1 1 — I I I I 1 1 1 



H 1 I I I I I I 



> 



V 



V 



J I I LJ_ 



0.1 



10 



PSD Break Time Scale (days) 



Fig. 12. — Plot of bolomctric luminosity and estimated black hole mass plotted against characterstic variability time scale (reciprocal of 
the PSD break frequency) from the singly-broken power law model fits (for NGC 3783, the high-frequency break from the doubly-broken 
model fit is used; lower limits to time scale are shown for Fairall 9 and NGC 5548). The time scale for MCG-6-30-15 is from UMP02. Masses 
are reverberation mapping estimates from Kaspi et al. (2000) except for NGC 3516 (from Wanders & Home 1994) and Akn 564 (rough 
reverberation estimate by Collier et al. 2001). MCG— 6-30-15 does not have a reverberation mass estimate; see text for details. The dashed 
line denotes the linear mass-time scale relation T( days ) = Mbh/10 6 ' 5 Mq. 
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Fig. 13. — Plot of bolometric luminosity and estimated black hole mass plotted against the amplitude of the best-fitting underlying PSD 
at the break frequency. The model fit and amplitude for MCG-6-30-15 is from UMP02. 



